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We present an update on Numerical Stochastic Perturbation Theory projects for Lattice QCD, 
which are by now run on apeNEXT. As a first issue, we discuss a strategy to tackle finite size 
effects which can be quite sizeable in the computation of logarithmically divergent renormal- 
ization constants. Our first high loop determination of quark bilinears for Wilson fermions was 
limited to finite constants and finite ratios. A precise determination of Zp and (and hence of 
Z,„) now becomes possible. We also give an account of computations for actions different from 
the standard regularization we have taken into account so far (Wilson gauge action and Wilson 
fermions). In particular, we present the status of computations for the Lattice QCD regularization 
defined by tree level Symanzik improved gauge action and Wilson fermions, which became quite 
popular in recent times. We also take the chance to discuss the related topic of the computation of 
the gluon and ghost propagators (which we undertook in collaboration with another group). This 
is relevant in order to better understand non-perturbative computations of propagators aiming at 
qualitative/quantitative understanding of confinement. 
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1. Introduction 



Numerical Stochastic Perturbation Theory (NSPT) proved to be a viable tool to compute high 
loop renormalization constants (RC) for Lattice QCD In the following we address a couple 
of the current issues in this subject. One is a problem: in our previous work we pointed out that 
finite volume effects can be quite important for quantities having an anomalous dimension (i.e. 
displaying log-divergencies). A second one is an opportunity: taking into account different actions 
(both for gluons and for fermions) does not imply a real overhead (there are no Feynman rules 
to derive). In the following we will report on these items, adding as an interesting aside a brief 
account on the gluon and ghost propagators, which, as a matter of fact, had not yet been looked at 
by NSPT 



2. Computation of renormalization constants 

RI-MOM' is a convenient scheme to compute RC by NSPT. A useful example is the compu- 
tation of quark bilinears. First of all we compute them (in Landau gauge) between external quark 
states at fixed (off-shell) momentum p 



J dx{p\ \if{x)T^f{x) \p) = Gr{pa) 



This step is taken in NSPT much the same way as in non-perturbative simulations. One then 
amputates {S{pa) is the quark propagator) and projects out the tree-level structure 

rr(pa) = S^^{pa) Gr{pa) S^^{pa) Or{pa) = Tr[Po^ Tripa)) 
Renormalization conditions now read (Z^ is the quark field renormalization constant) 



(2.1) 



3. Dealing with anomalous dimensions 

In NSPT one would like to take anomalous dimensions for granted. For RI-MOM' this is 
actually the case [Q]. Let us keep this in mind while writing our master formula for the scalar 
current 

(i-V+...)(i+^^^V^+...)(i+^^^+...) , , = 1 

p p p P^=M^ 

We explicitly wrote both the constant and the logarithmic contributions to renormalization con- 
stants (at first order the only log comes from Zj since one-loop quark-field anomalous dimension 
is zero in Landau gauge), ol^' {p^ 
can solve the previous relation to 



is zero in Landau gauge), ol^' (p^) is what is actually numerically measured. At one-loop order we 



z^^^-zl'^ = oY>{p')-rriog{p'). (3.1) 
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Figure 1: On the left: computation of one loop renormalization constant for the scalar current. In the 
notation of Eq. (3.1), upper points are the unsubtracted oi^\p^), while lower (circled crosses) stand for 
the subtracted oi^\p^) — 7i^'log(p^). Analytic result (darker symbol) is missed. On the right: the same 
computation subtracting the tamed-log contribution referred to in the text. Analytic result is got. 



The quantity on l.h.s. is now finite and the only dependence on pa = p is an irrelevant 
one, which can be wiped out by extrapolating to zero by means of what we call an Hypercubic- 
symmetric Taylor expansion. Unfortunately this program results in a failure: see Fig. (1) (left). 

IR problems do not come as a surprise, due to finite size effects. In Fig. (2) we display results 
on sizes 32^ and 16'*. Results perfectly agree for the ratio of Zg and Zp (which is a finite quantity). 
On the other hand, the smaller the lattice, the bigger the IR effects one gets for log-divergent 
quantities. These are one loop results, but the picture stays much the same at higher loops. 
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Figure 2: Computations of 0^p\p^) (the equivalent of Eq. (3.1) for the pseudoscalar current) (top) and 
(9l-''(/3^) (bottom) on 32^ (circles) and 16'* (diamonds). In the middle the ratio ^^p^^y 
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Circumventing the problem is not so hard at one loop. Define L = Na. Let us now write 
down for the quantity at hand the momentum sum I{p,a,L) of conventional Lattice Perturbation 
Theory, making use of the same regularization of zero modes that we use in NSPT (zero momentum 
removed from the sum). In the same spirit of ^ let us now split it as 

I{p,a,L) = I{0,a,L) + {I{p,a,L) - I{0,a,L)) = I{0,a,L) +J{p,a,L). (3.2) 

The divergence is logarithmic so that by subtracting 7(0, a, L) we make J{p,a,L) UV finite. IR 
divergences will pop up and cancel with those in 1(0, a, L): 

I{0,a,L) = Ci + ylog{a/L)+H{a/L) (3.3) 
J{p,a,L) = C2 + ylog{pL) + G{pa,a/L,pL) . 

Now we look for pL = pN effects. These should be looked for in G{pa,a/L,pL) — G{pL), which 
we computed from the formal continuum limit of our sum J{p,a,L): this means a — > with L = Na 
fixed (as we pointed out, J{p,a,L) is UV finite). Again, we need the same ad hoc regularization 
of zero modes. We call this contribution a tamed-log. It is supposed to resemble the expected log, 
but with pL = pN effects on top of it: it indeed approaches a log for >> 1. By subtracting this 
tamed-log we got the right one loop result. 

4. Dealing with anomalous dimensions: a new strategy. 

The above method has got obvious drawbacks. First of all, one has to go back to diagramatic 
computations we would have liked to get rid of by means of NSPT. While one is happy enough 
with the one loop picture, at higher loops the situation is less clear and it is for sure much more 
cumbersome. Moreover, one is not actually making use of the computations on different lattice 
sizes, but has to revert instead to a continuum computation. In the following we present a strategy 
which is currently under investigation (we do not yet display results). This method is a very close 
relative of the one which is employed in an NSPT 3-d application with a mass in play [Q]. 

It can be better understood having in mind figure 2. An obvious way of looking at it is the 
following: on 32^ and 16^ one gets results at the same physical momentum p from points affected 
by different pL = pN effects. Moreover, for pL large enough one can get results substatially free 
of finite volume effects. Having measurements on at least two different sizes (A^i, N2, Ni > N2), 
this suggests to proceed as follows: 

• Go to a momentum p high enough: the measurement f{p,pN\ » 1) ^ f{p,°°) on the big 
lattice will be substantially free of finite size effects. 

• Read the deviation 5{pN2) = f{p,pN2) —f{p,°°) ~ f{p,pN2) — f{p,pN[) {i.e. the deviation 
from the result you got on the little lattice at the same value of p). 

• This deviation defines the pL = pN effect you have to correct for at the point having the same 
value of pL = pN. 

The strategy is promising and first results are encouraging. 



4 



High loop renormalization constants by NSPT: a status report 



Francesco Di Renzo 



1.9r 



-0.48 r 



-0.52 



E-1.7- 







1.6 



0.2 



0.4 0.6 
(pa)' 



-0.56 



-0.6 L 



0.2 



0.4 0.6 

(pa)' 



Figure 3: Computations of the Wilson fermions one loop critical mass (left) and self-energy (right) for 
tree-level Symanzik improved gauge action. For the self-energy points come in "families" corresponding to 
different violations of Lorents symmetry, depending on the lenght of momentum along the gamma matrix 
one traces with. Results from the measurement on a single configuration. See text for an important caveat. 

5. NSPT computations for different actions 

Changing the action can be quite cumbersome for Lattice Perturbation Theory: deriving Feyn- 
man rules for propagators and (both relevant and irrelevant) vertices can require hard work (this is 
usually a computer-aided task). Life is by far easier in NSPT. We are working on different Lattice 
QCD regularizations resulting from taking different gauge and fermionic actions. For the gauge 
action we consider plain Wilson, Iwasaki or (tree-level) Symanzik. For fermionic action we take 
plain Wilson or Clover. From the computational point of view, notice that all the combinations are 
fairly well implemented on apeNEXT (as compared to APEmille) because of the larger amount of 
available memory. 

Fig. (3) displays one loop benchmark computations of critical mass and fermion self-energy 
for the tree-level Symanzik improved gauge action with plain Wilson fermions. At the time of 
the conference we were still in a bit embarassing situation: from Fig. (3) one can inspect results 
were (slightly) off the analytic results available in the literature [||, |6|] (a wrong normalization is 
now being fixed). Notice with this respect that a one loop benchmark is a very safe one for NSPT: 
since everything is implemented order by order, once you get first order right you are standing on 
a quite safe ground. Apart from the actual results, notice that plots are measurements from a single 
configuration. Fermionic measurements are remarkably stable (curves are smooth enough): keep 
this in mind now that we move to the Gluon Propagator. 



6. Gluon (and Ghost) Propagator 

Gluon and ghost propagators are the subject of a new collaboration^ . The framework for this 
work is the qualitative/quantitative investigation of confinement via Schwinger-Dyson equations 
for the fundamental degrees of freedom of the theory. 



'F. Di Renzo, C. Torrero, M. Ilgenfritz, M. MuUer-Preussker, H. Perlt, A. Schiller. 
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Figure 4: Tree level (left) and one loop (right) computation of the gluon propagator. Notation is: Ylnvik) = 

^{k) p"^ — ■ The analytically known constants are recovered in the continuum limit (at one loop the log 

dictated by anomalous dimension has to be subtracted). Also in this case data arrange themselves in families 
according to different patterns of violation of Lorentz symmetry. 

The gluon propagator in NSPT is a really straightforward computation. One first fixes the gauge: 
once again, Landau is the obvious choice (FFT acceleration is a powerful tool at hand). One then 
goes to momentum space, takes the product of fields and traces. The ghost propagator is measured 
much the same way as the fermion propagator, by inverting order by order the Faddeev-Popov 
matrix on a (momentum space) source. The technology for the ghost was derived in the context of 
another application [^. 

In figure 4 we display tree-level and one loop benchmark computations of the gluon propa- 
gator. Notice that this time we present also a tree level computation as a benchmark. The gluon 
propagator is with this respect quite different from the fermion (or ghost) propagator: there is no 
inversion on a source, but simply the computation of a correlator (in momentum space). While this 
is straightforward, one now gets much more noise. Even tree level comes from an actual measure- 
ment: first fluctuations around the vacuum are the Lie algebraic free fields whose correlatator gives 
the free Feynman propagator. It is funny to compare the previous plots (which come out of O(IOO) 
configurations) with the (one configuration) measurements of Fig. (3). Still, although noisy, this is 
a doable computation. Both for the gluon and for the ghost, we provide measurements from fairly 
big lattices (16^ and 32^). These measurements are taken on the same configurations stored to mea- 
sure the QCD fermionic quantities. Our German collaborators [||] provide instead huge statistics 
measurements from smaller lattices. 

7. Conclusions 

We have presented a new method to correct for finite volume effects in NSPT renormalization 
constants computations. The method is much less cumbersome than the one we applied to one loop 
computations. 
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We are computing renormalization constants for different regularizations of Lattice QCD: 
computations for Iwasaki and tree level Symanzik improved gauge actions are on their way. 

We also reported on the computation of the gluon (and ghost) propagator. While the gluon 
propagator is the most straightforward to measure, it is not at all the best one with respect to the 
statistical fluctuations. 
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